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Abstract 

While many dynamical systems of mechanical origin, in particular 
billiards, are strongly chaotic - enjoy exponential mixing, the rates of 
mixing in many other models are slow (algebraic, or polynomial). The 
dynamics in the latter are intermittent between regular and chaotic, 
which makes them particularly interesting in physical studies. How- 
ever, mathematical methods for the analysis of systems with slow mix- 
ing rates were developed just recently and are still difficult to apply to 
realistic models. Here we reduce those methods to a practical scheme 
that allows us to obtain a nearly optimal bound on mixing rates. We 
demonstrate how the method works by applying it to several classes 
of chaotic billiards with slow mixing as well as discuss a few examples 
where the method, in its present form, fails. 



1 Introduction 

A billiard is a mechanical system in which a point particle moves in a compact 
container Q and bounces off its boundary dQ. It preserves a uniform measure 
on its phase space, and the corresponding collision map (generated by the 
collisions of the particle with dQ, see below) preserves a natural (and often 
unique) absolutely continuous measure on the collision space. The dynamical 
behavior of a billiard is determined by the shape of the boundary dQ, and it 
may vary greatly from completely regular (integrable) to strongly chaotic. 
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In this paper we only consider planar billiards, where Q C M^. The 
dynamics in simple containers (circles, ellipses, rectangles) are completely 
integrable. The first class of chaotic billiards was introduced by Ya. Sinai 
in 1970 [25; he proved that if dQ is convex inward and there are no cusps 
on the boundary, then the dynamics is hyperbolic (has no zero Lyapunov 
exponents), ergodic, mixing and K-mixing. He called such systems dispersing 
billiards, now they are often called Sinai billiards. Gallavotti and Ornstein 
proved in 1976 that Sinai billiards are Bernoulli systems. 

Many other classes of planar chaotic billiards have been found by Buni- 
movich OH] in the 1970s, and by Wojtkowski [22], Markarian ^H], Donnay 
|13j . and again Bunimovich [H] in the 1980s and early 1990s. All of them are 
proven to be hyperbolic, and some - ergodic, mixing, and Bernoulli. 

Multidimensional chaotic billiards are known as well, they include two 
classical models of mathematical physics - periodic Lorentz gases and hard 
ball gases, for which strong ergodic properties have been established in a 
series of fundamental works within the last 20 years, we refer the reader to 
a recent collection of surveys |15j . 

However, ergodic and mixing systems (even Bernoulli systems) may have 
quite different statistical properties depending on the rate mixing (the rate 
of the decay of correlations), whose precise definition is given below. On 
the one hand, the strongest chaotic systems - Anosov diffeomorphisms and 
expanding interval maps - have exponential mixing rates, see a survey ^2]- 
This fact implies the central limit theorem, the convergence to a Brownian 
motion in a proper space-time limit, and many other useful approximations 
by stochastic processes that play crucial roles in statistical mechanics. 

On the other hand, many hyperbolic, ergodic and Bernoulli systems have 
slow (polynomial) mixing rates, which cause weak statistical properties. Even 
the central limit theorem may fail, see again a survey ^2]- Such systems are, 
in a sense, intermittent, they exemplify a delicate transition from regular 
behavior to chaos. For this reason they have attracted considerable interest 
in mathematical physics community during the past 20 years. 

The rates of mixing in chaotic billiards is rather difficult to establish, 
though, because the dynamics has singularities, which aggravate the analysis 
and make standard approaches (based on Markov partitions and transfer 
operators) inapplicable. For planar Sinai billiards, L.-S. Young [21| developed 
in 1998 a novel method (now known as Young's tower construction) to prove 
exponential (fast) mixing rates. Young applied it to Sinai billiards under 
two restrictions - no corner points on the boundary and finite horizon. Her 
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method was later extended to all planar Sinai billiards jHj and to more general 
billiard-like Hamiltonian systems fOJ P . 

There are no rigorous results on the decay of correlations for multidi- 
mensional chaotic billiards yet, since even the best methods cannot handle 
a recently discovered phenomenon characteristic for billiards in dimension 3 
and higher - the blow-up of the curvature of singularity manifolds [2] . 

For planar billiards with slow (nonexponential) mixing rates, very little 
is known. They turned out to be even harder to treat than Sinai billiards. 
Most published accounts are based on numerical experiments and heuristic 
analysis, which suggest that Sinai billiards with cusps on the boundary [T7] 
and Bunimovich billiards [22] have polynomial mixing rates. Slow mixing in 
these models appears to be induced by tiny places in the phase space, where 
the motion is nearly regular, and where the moving particle can be caught 
and trapped for arbitrarily long times. 

A general approach to the studies of abstract hyperbolic systems with 
slow mixing rates was developed by Young in 1998 |21], but its application to 
billiards required substantial extra effort, and it took five more years. Only in 
2004 R. Markarian ^H] used Young's method to establish polynomial mixing 
rates for one class of chaotic billiards - Bunimovich stadia. He showed that 
the correlations for the corresponding collision map decayed as 0(^n~^ In^ n) . 

In this paper we generalize and simplify the method due to Young and 
Markarian essentially reducing it to one key estimate that needs to be verified 
for each chaotic billiard table. This estimate, see (j6.2|) in Section ^ has 
explicit geometric meaning and can be checked by rather straightforward 
(but sometimes lengthy) computations, almost in an algorithm-like manner. 
We also remark that the method is not restricted to billiards, it is designed 
for general nonuniformly hyperbolic systems with Sinai-Ruelle-Bowen (SRB) 
measures. 

Then we apply the above method to several classes of billiards. These 
include semi-dispersing billiards in rectangles with internal scatterers, Buni- 
movich flower-like regions, and "skewed" stadia ("drive-belt" tables). In 
each case we prove that the correlations decay as O (^n~°' In^^'^ n) for a cer- 
tain a > 0, which depends on the character of "traps" in the dynamics, and 
is different for different tables. We also show that for some billiard tables 
our key estimate fails, hence the correlation analysis requires further work. 
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2 Statement of results 



Here we state our estimates on the decay of correlations for several classes of 
billiard systems. The general method for proving these estimates is described 
in the next two sections. 

First we recall standard definitions [HI El El Ej- A billiard is a dynamical 
system where a point moves freely at unit speed in a domain Q [the table) 
and reflects off its boundary dQ {the wait) by the rule "the angle of inci- 
dence equals the angle of reflection". We assume that Q C and dQ is a 
finite union of curves (arcs). The phase space of this system is a three 
dimensional manifold Q x S^. The dynamics preserves a uniform measure 
on g X S\ 

Let Ai = dQ X [— 7r/2,7r/2] be the standard cross-section of the billiard 
dynamics, we call Ai the collision space. Canonical coordinates on A4 are 
r and ip, where r is the arc length parameter on dQ and (f G [— 7r/2,7r/2] is 
the angle of reflection, see Fig. [T] 




(a) (b) 

Fig. 1: Orientation of r and ip 

The first return map : Ai Ai is called the collision map or the 
billiard map, it preserves smooth measure dfi = cos ipdr dip on Ai. 
Let f,g& L'jj^{Ai) be two functions. Correlations are defined by 

(2.1) Cn{f,g,J',fi)= [ {foJ^-)gdfi- [ fdfif gdfi 

Jm Jm Jm 

It is well known that JF : ^ 7VI is mixing if and only if 

(2.2) limC„(/,^?,^,/i) = \/f,geLl{M) 

The rate of mixing of JF is characterized by the speed of convergence in ()2.2p 
for smooth enough functions / and g. We will always assume that / and g 
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are Holder continuous or piecewise Holder continuous with singularities that 
coincide with those of the map for some k. For example, the free path 
between successive reflections is one such function. 
We say that correlations decay exponentially if 

|C„(/,fi', < const ■ e"''" 

for some c > and polynomially if 

\Cn{f,g,J^,^)\ < const -n"" 

for some a > 0. Here the constant factor depends on / and g. 
Next we state our results. 

Let 7^ C be a rectangle and . . . , C int 7?. open strictly convex 
subdomains with smooth or piecewise smooth boundaries whose curvature is 
bounded away from zero and such that Bi fl ,Bj = for i ^ see Fig. |2] (a). 
A billiard 'uiQ = TZ \ UiBi is said to be semi- dispersing since its boundary is 
partially dispersing (convex) and partially neutral (fiat); the flat part is dTZ. 

Theorem 1. For the above semi- dispersing billiard tables, the correlations 
\2. 1]] for the billiard map T : M. ^ M. and piecewise Holder continuous 
functions f,g on Ai decay as |C„(/, (7, JF, /i)| < const ■ (lnn)^/n. 




(a) (b) (c) 

Fig. 2: Slow mixing billiard tables 

Next, let Q C be a domain with a piecewise smooth boundary 

= r = Fi U ■ ■ • u r, 

such that each smooth component F.j C dQ is either convex inward (dispers- 
ing) or convex outward (focusing). Assume that every focusing component 
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Fj is an arc of a circle such that there are no points of dQ on that circle or 
inside it, other than the arc Fj itself, see Fig. |21 (b). Billiards of this type 
were introduced by Bunimovich j3] who established their hyperbolicity and 
ergodicity. 

We make two additional assumptions: (i) if two dispersing components 
Fj, Fj have a common endpoint, then they are transversal to each other at 
that point (no cusps!); (ii) every focusing arc Fj C dQ is shorter than half 
a circle. We also assume that the billiard table is generic to avoid certain 
technical complications, see Section |H| 

Theorem 2. For the above Bunimovich-type billiard tables, the correlations 
\2. 1]) for the billiard map T : M. ^ M. and piecewise Holder continuous 
functions f,g on Ai decay as |C„(/, (7, JF, < const ■ (InnY/n^. 

In fact, the boundary dQ may have flat components as long as there is an 
upper bound on the number of consecutive reflections off the fiat components. 
For example, this is the case when dQ has a single fiat component or exactly 
two nonparallel fiat components. Then the above theorem remains valid. 
However, if the billiard particle can experience arbitrarily many consecutive 
collisions with fiat boundaries, then, generally speaking, even the ergodicity 
cannot be guaranteed, let alone mixing or decay of correlations. 

Lastly, we consider a special stadium. It is a convex domain Q 

bounded by two circular arcs and two straight lines tangent to the arcs at 
their common endpoints, see Fig. |21(c). We distinguish between a "straight" 
stadium, whose fiat sides are parallel, see Fig. |21(c), and a "skewed" stadium, 
whose fiat sides are not parallel, we call them drive-belt tables, see Fig.lTUlfa). 
Notice that drive-belt tables contain an arc longer than a half circle, so they 
do not satisfy the assumption (ii) of the previous theorem (or the assumptions 
made by Markarian in jTHj). Again billiards of this type were introduced by 
Bunimovich |1| who established their hyperbolicity and ergodicity. 

Theorem 3. For both types of stadia, the correlations \2.1\] for the billiard 
map T : M. ^ M. and piecewise Holder continuous functions f,g on M. 
decay as |C„(/, (7, JF, /i)| < const ■ (lnn)^/n. 

The same bound on correlations for straight stadia is already obtained 
by Markarian ^H]; but we include it here for the sake of completeness. 



6 



3 Correlation analysis 



Here we present a general method for estimating correlations in nonuniformly 
hyperbolic dynamical systems. It is based on Young's recent results [211 
and their extensions by Markarian UHj and one of us [Oj. 

Let JF : ^ be a nonuniformly hyperbolic map acting on a Rieman- 
nian manifold A4 with or without boundary. We assume that J-' preserves 
an ergodic Sinai-Ruelle-Bowen (SRB) measure fi, see E] for definitions 
and basic facts. 

Young [211 123 found sufficient conditions under which correlations for the 
map decay exponentially. We only sketch Young's method here, skipping 
many technical details, in order to focus on the main condition, see (j3.1|) 
below. 

The key element of Young's construction is a "horseshoe" Aq (a set with 
a hyperbolic structure, often called a rectangle). By iterating points x G Ag 
under the map J-" until they make proper returns to Aq, see definitions in [21], 
Young constructs a tower. A, in which the rectangle Aq constitutes the first 
level (the base). The induced map J^a on the tower A moves every point one 
level up until it hits the ceiling upon which it falls onto the base Aq. Now, 
for X G A, let _R(x; JF, Aq) = mm{k > 1 : JF|(x) G Aq} denote the "return 
time" of the point x to the base Aq. The tower has infinitely many levels, 
hence R is unbounded. Young proves that if the probability of long returns 
is exponentially small, then correlations decay exponentially fast. Precisely, 
if 

(3.1) G A : R{x; Aq) > n) < const ■ ^" Vn > 1 

where 6' < 1 is a constant, then |Cn(/, (?, T \ < const ■ e~'^" for some c > 0. 

Therefore, if one wants to derive an exponential bound on correlations 
for a particular map, one needs to construct a "horseshoe" Aq and verify 
the tail bound p.ip . The latter involves all the iterates of the map JF. To 
simplify this verification one of us (NC) proposed P sufficient conditions, 
which involve just one iterate of the map JF (or one properly selected power 
!F"^, see the next section), that imply the tail bound ()3.1|) . Those conditions 
do not require a construction of a horseshoe Aq. We present those conditions 
fully in the next section. 

Young later extended the results of [21] to cover system with slow mixing 
rates. She proved [22] that if 

(3.2) i2{xe A: R{x; Aq) > n) < const ■ n'^ \/n > 1 
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where a > is a constant, then correlations decay polynomially: 



(3.3) 



Cn{f,g,^,fJ')\ < const -n 



—a 



Again, a direct verification of the polynomial tail bound ()3.2p in specific 
systems involves all the iterations of the map JF, and thus might be difficult. 
There are no known simplifications, like the one mentioned above, that would 
reduce this problem to just one iterate of the map either. However, there is 
a roundabout way proposed by Markarian ^H] which simplifies the analysis, 
even though it leads to a slightly less than optimal bound on correlations. 
We describe Markarian's approach next. 

Suppose one can localize places on the manifold Ai where the dynamics 
fails to be strongly hyperbolic, and find a subset M C on which has a 
strong hyperbolic behavior. This means, precisely, that the first return map 
F : M ^ M satisfies the conditions of Young's paper |21]; in particular, 
there exists a horseshoe Aq C M for which an exponential tail bound (|3.1|) 
holds for the return times under the map F. Now, of course, the verification 
of the bound (j3.1|) for the map F can be done via the simplified conditions 



Next, consider the return times to M under the original map i.e. 



for X G A^. Suppose they satisfy a polynomial tail bound 

(3.5) ^i{xe M: R{x; M) > n) < const ■ n"" Vn > 1 
where a > is a constant. Equivalently, we may suppose that 

(3.6) /i(x e M : R{x; M) > n) < const ■ n~"~^ \/n > 1 

The equivalence of (j3.5p and (|3.6|) is a simple fact, we omit its proof. Tail 
bounds for the return times R{x; JF, M) to M are much easier to obtain, 
in many systems, than those for the return times R{x] JF, Aq) to Aq. The 
following theorem was essentially proved in , but we provide a proof here 
for completeness: 

Theorem 4. Let : M. ^ M. he a nonuniformly hyperbolic map. Suppose 
M d M. is a subset such that the first return map F : M M satisfies the 
tail bound \3. 1]) for the return times R{x; F, Aq) to a rectangle Aq C M. If the 



of 0. 



(3.4) 



R{x]J^, M) 



min{r > 1 : JF'^(x) G M} 
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return times R{x; JF, M) satisfy the polynomial bound i3. 5)) or, equivalently, 
larm . then 

(3.7) |C„,(/,^,J^,/i)| < const ■ {InnY^^n-'' 
Proof. For every n > 1 and x & M. denote 

r{x;n,M) = #{1 < i < n : r{x) G M} 

and 

An = {x e M: R{x; Ao) > n}, 
-Bn,6 = {xGA^: r(x; n, M) > 61nn}, 

where 6 > is a constant to be chosen shortly. By ()3.H) . 

/i(v4„n5„,5) < const -n^^'"^". 

Choosing h large enough makes this bound less than const-n~". 

To bound [li^A^ \ Bn^b) we note that points x & An\ B^^ return to M 
at most 6 In 77, times during the first n iterates of JF. In other words, there 
are < 61nn time intervals between successive returns to M, and hence the 
longest such interval, we call it /, has length > n/(61nn). Applying the 
bound ()3.6|) to the interval / gives 

/i(A„ \ Bn,b) < const ■ n (lnn)"+V^"^^ 

(the extra factor of n must be included because the interval I may appear 
anywhere within the longer interval [l,n], and the measure n is invariant). 
In terms of Young's tower A, we obtain 

(3.8) /i(x e A : R{x; Aq) > n) < const ■ (lnn)'^+^ n"'' Vn > 1 

This tail bound differs from Young's (|3.2p by the extra factor (Inn)"^^. How- 
ever, as it was explained by Markarian ^2]! the same argument that Young 
used to derive ()3.3p from ()3.2|) now gives us ()3.7|) based on ()3.8|) . This com- 
pletes the proof of the theorem. □ 
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4 Conditions for exponential mixing 



Here we list sufficient conditions on a 2D nonuniformly hyperbolic map 
F: M — >■ M with a mixing SRB measure fi under which its correlations de- 
cay exponentially. These are a 2D version of more general conditions stated 
in [Hj. We also provide comments that will help us apply these conditions to 
chaotic billiards. 

We assume that M is an open domain in a two-dimensional C°° compact 
Riemannian manifold Ai with or without boundary. For any smooth curve 
W C ^A we denote by the Lebesgue measure on W (induced by the 
Euclidean metric). For brevity, = uwi^) will denote the length of W. 

4.1 Smoothness. The map F is a diffeomorphism of M\S onto F{M\S), 
where 5 is a closed set of zero Lebesgue measure. Usually, S is the set of 
points at which F either is not defined or is singular (discontinuous or not 
differentiable) . 

Remark. The collision map J^: M. ^ M. for a billiard table with a piece- 
wise smooth boundary is piecewise C""^^ smooth The singularities of 
J-' make a closed set S that is a finite or countable union of smooth compact 
curves in the collision space M.. On those curves, one-sided derivatives of JF 
are often infinite. When we construct a ffist return map F: M — > M on a 
subdomain M C it will be always clear that F has similar properties. 

We denote by Sm = S U J^~^{S) U • ■ ■ U .F-™+^(5) the singularity set for 
the map JF*". Similarly, S^m denotes the singularity set for the map JF"™-. 

4.2 Hyperbolicity. There exist two families of cones (unstable) and 
(stable) in the tangent spaces TxM, x G M, such that DF{C^) C Cp^ and 
DF{Cl) D Cp^ whenever DF exists, and 

\DF{v)\>K\v\ 'iveCl and \DF-\v)\> \\v\ 'iv e 

with some constant A > 1. Here | ■ | is some norm equivalent (though not 
necessarily uniformly equivalent) to the Euclidean norm. These families of 
cones are continuous on M and the angle between and C| is bounded 
away from zero. 

For any F-invariant measure //', almost every point a; G M has one pos- 
itive and one negative Lyapunov exponent. Also, almost every point x has 
one-dimensional local unstable and stable manifolds, which we denote by 
W^i^x) and IV (x), respectively. 
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Remark. The existence of Lyapunov exponents usually follows from the 
Oseledec theorem, and the existence of stable and unstable manifolds - from 
the Katok-Strelcyn theorem, see fH]- Both theorems require certain mild 
technical conditions, which have been verified for virtually all planar billiards 
in jni- The hyperbolicity has been proven for all the classes of billiards 
studied in this paper. It is also known that the tangent vectors to the curves 
in S lie in stable cones, and those of the curves in - in unstable cones. 

4.3 SRB measure. The map F preserves an ergodic measure /x whose 
conditional distributions on unstable manifolds are absolutely continuous. 
Such a measure is called a Sinai-Ruelle-Bowen (SRB) measure. We also 
assume that /i is mixing. 

Remark. For the collision map of a chaotic billiard system, the natural 
smooth invariant measure /i is an SRB measure [21]. Its ergodicity and 
mixing have been proven for all classes of billiards studied in this paper, see 
references in Section |21 

4.4 Distortion bounds. Let A(a;) denote the factor of expansion on the 
unstable manifold W^lx) at the point x. If x,y belong to one unstable 
manifold such that F"- is defined and smooth on W^, then 

(4.1) ^ogll^^^<^dist{F-x,F-y)) 

where ■?/'(■) is some function, independent of W, such that ip{s) — > as s ^ 0. 

Remark. Since the derivatives of the collision map turn infinite on 
some singularity curves S, the expansion of unstable manifolds terminating 
on those curves is highly nonuniform. To enforce the required distortion 
bound, one follows a standard procedure introduced in [Zj, see also PHE]- 
Let M.+ denote the part of the collision space corresponding to dispersing 
walls Ti C dQ, i.e. 

^A+ = {{r,{p) : r G Fj, Fj is dispersing}. 

One divides Ai+ into countably many sections (called homogeneity strips) 
defined by 

Hk = {{r,ip) e M+: tt/2 - k-"^ < if < n/2 - {k + 1)-^} 
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and 



= {{r,^)eM+: - 7i/2 + {k + 1)"^ <^< -n/2 + k~^} 
for all k > ko and 

(4.2) Ho = {{r,^)eM+: - n/2 + kf < ^ < 7r/2 - kf}, 

here /cq > 1 is a fixed (and usually large) constant. Then, a stable (unstable) 
manifold W is said to be homogeneous if its image J^^iW) lies either in one 
homogeneity strip of or in \ for every n > (resp., n < 0). It 
is shown in |7] that a.e. point x E Ai has homogeneous stable and unstable 
manifolds passing through it. The distortion bounds for homogeneous unsta- 
ble manifolds were proved in [7[ E] . It is also easy to check that the distortion 
bounds for the collision map imply similar bounds, with the same function 
ip{s), for the induced map F : M — > M on any subdomain M C M.. 

iFiom. now on, we will only deal with homogeneous manifolds without 
even mentioning this explicitly. This can be guaranteed by redefining the 
collision space, as it is done in we remove from M.^ the boundaries dH^, 
thus making a countable disjoint union of the open homogeneity strips 
iffc's. (There is no need to subdivide or otherwise redefine ^A \ A^+, see P) 
Accordingly, the images (preimages) of dHk need to be added to the set S-i 
(resp. S). We call them new singularities, and refer to the original sets S 
and 5_i as old singularities. Now the stable and unstable manifolds for the 
map JF on the (thus redefined) collision space M. will be always homogeneous 
and the distortion bounds will always hold. 



4.5 Bounded curvature. The curvature of unstable manifolds is uniformly 

bounded by a constant i? > 0. 

4.6 Absolute continuity. If Wi, W2 are two small unstable manifolds close 
to each other, then the holonomy map h : Wi W2 (defined by sliding 
along stable manifolds) is absolutely continuous with respect to the Lebesgue 
measures and h'w2y ^.nd its Jacobian is bounded, i.e. 

(4.3) l/C < < C 

with some C > 0, where Wl C Wi is the set of points where h is defined. 
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Remark. The properties 4.5 and 4.6 have been verified for planar chaotic 
bilhards considered in this paper [7| |H] and even for more general billiard-like 
Hamiltonian systems [TU] . 

Before we state the last (and most important) condition, we need to 
introduce some notation. Let 60 > 0- We call an unstable manifold W a 60- 
LUM (local unstable manifold) if \W\ < 6q. Let V G W he an open subset, 
i.e. a finite or countable union of open sumbintervals of W. For x G V denote 
by V{x) the subinterval of V containing the point x. Let n > 0. We call an 
open subset V CW a. {6q, n)-subset if the map F" is defined and smooth on 
V and |F"y(x)| < 60 for every x & V. Note that F"V^ is then a union of 
(5o-LUM's. Define a function ry„ on V by 

(4.4) ry,„(a;) = dist(F"x,(9F"\/(x)), 

which is the distance from F'^{x) to the nearest endpoint of the curve F'^V{x). 
In particular, twq^x) = dist{x,dW). 

4.7 One-step growth of unstable manifolds. There is a small 6q > and 
constants G (0, 1) and j3o, D, K,a > with the following property. For any 
sufficiently small 6 > and any (5o-LUM W denote by = Ug{W) C W the 
5-neighborhood of the set W (1 S within W. Then there is an open {Sq, 1)- 
subset = Vs\W) CW\U} such that uw{W \ (f/j U 1//)) = and V5 > 

(4.5) ^w{rvi,i < ^) < 2aoe + e/9o<5(r^|"W^| 



(4.6) ^w{ruj,o<£) <D6~^e 
and 

(4.7) i^w{Ul)<D6'^ 

We comment on these conditions after stating the theorem proved in jH] : 

Theorem 5 (jHl)- Let F satisfy the assumptions 4-1^4-'^- Then (a) there is 
a horseshoe Aq C M such that the return times R{x] F, Aq) satisfy the tail 
bound i3. 1]) . Thus, (h) the map F : M —>■ M enjoys exponential decay of 
correlations. 
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Remark. It is shown in jH], Proposition 10.1, that if the conclusions (a) 
and (b) of this theorem hold for some power F™, m > 2, of the map F, then 
they hold for F itself. Therefore, it will be enough to verify the condition 
4.7 for the map F"^ with some m > 2 (of course, S should then be replaced 
by Sm, by ry^m, etc.). This option will be important when the condition 
(14. 5|) fails for the map F but holds for some power F"^, see below. In this 
case we do not have to worry about (j4.6|) and (jUTj), since they will hold for 
all F™, m > 1, see the next section. 

5 Simplified conditions for exponential mix- 
ing 

Here we simplify the conditions ()4.5p - ()4.7|) and reduce them to one inequality 
that needs to be checked for every particular class of chaotic billiards. 

First, the condition (j4.6|) always holds for 2D maps, where dimH^ = 1, 
as it was shown in jU], and the reason is simple: each connected component 
of has length > 2S, thus there are no more than 60/ (26) of them, hence 
^w^i^ulfi < ^) < SoS^^e. Thus, fl4.6|) will hold for all our applications. 

The condition ()4.7|) holds under very mild assumptions on the structure 
of the set S, which will be always satisfied in our applications: 

Assumption (Structure of the singularity set). Assume that for any 
unstable curve W G M (i.e., a curve whose tangent vectors lie in unstable 
cones) the set W (1 S is finite or countable and has at most K accumulation 
points on W, where K > 1 is a constant. Let Xoo be one of them and 
denote the monotonic sequence of points WnS converging to Xoo- We assume 
that 



for some constants C, > 0, i.e. the convergence of Xn to x^o is faster than 
some power function of n. 

Lemma 6. In the notation of 4-7, our assumption i5.1]} implies 



(5.1) 



dist(x„,Xoo) <C/n'^ 



(5.2) 



which in turn implies d^. 7| ). 
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1 ]_ 

Proof. It is obviously enough to consider the case K = 1. Let n^, = C i+<^5 i+f* 
and note that UjiW) is contained in the union of the interval (xoca^n.) and 
intervals of length 25 centered on Xj, 1 < i < n,,,, thus 

vw{Ul{W)) <C/nt + 2nJ 

which easily implies ()5.2|) . □ 

Lemma 7. Under our assumption \5.1\) . the hound d^. 7| j holds for any power 
(with D and a depending on m, of course). 



Proof. We use induction on m. Let W be an unstable manifold and V^iW) C 
W denote the 5-neighborhood of the set W fl Sm within W . Assume that 
(gm) holds, i.e. 

(5.3) vw{UT{W)) <D5'' 

for some constants -D, o" > and any 5 > 0. Set 

then i^w{Usl{W)) < D6l by ()5.3|) . Let Ji, . . . , Jat be the connected compo- 
nents oiW \Sm whose length is > 5*. Note that < < 6q/6^. For 
each i = 1, . . . ,N the set Ji = F"^{Ii) is an unstable manifold. By the as- 
sumption (|5.ip . the set JjfliS is at most countable and has < K accumulation 
points which satisfy 1)5.11) . Since is smooth and contracting on each J,, 
the set li n Sm+i also is at most countable and has < K accumulation points 
which satisfy 1)5.11) . Now by the same argument as in the proof of Lemma El 

for each i. Observe that 
thus 

iyw{Ur^\W)) < D6: + const ■ 6-^6^^ 
< const ■ 5(1+^11+^, 
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which proves (|4.7p for the map F"^~^^. □ 

Remark S. We need to extend Lemma IHl to a shghtly more general type of 
singularities. Suppose Si and S2 are two sets such that each satisfies the 
above assumption ()5.1|) . and S = SiU F~^{S2)- Then the argument used in 
the proof of Lemma [7| shows that S satisfies ()4.7p . Furthermore, in this case 
Sm satisfies ()4.7|) due to the same argument. 

We now turn to the most important condition (j4.5p . Let Wi denote the 
connected components of W \Si and 

A,- = min A(x) 

the minimal local expansion factor of the map F on Wi. We note that due 
to the distortion bounds 

(5.4) Vx, yew, e-^^'^^ < ^ < e^^''\ 

My) 

and so due to the smallness of 6q the expansion factor is in fact almost 
constant on each Wi. 

Lemma 8. The condition (|^.5| ) is equivalent to 

(5.5) «! : = liminf sup y ^A'[^ < 1, 

■50^0 W:\W\<So~' 

where the supremum is taken over unstable manifolds W and Aj, i > 1, 
denote the minimal local expansion factors of the connected components of 
W \ S under F. 

We call (j5.5j) a one-step expansion estimate for F. 

Proof. We prove that (j5.5|) implies (|4.5|) . the converse implication will be self- 
evident in the end (the converse will not be used in practical applications). 

Let Wl = Wi\U^ (W) be the open subinterval of Wi obtained by removing 
from Wi the 5-neighborhood of its endpoints (of course, if \Wi\ < 26, then 
Wl = 0), and let = UiWl. It is easy to see that We > 

(5.6) Jyw{rw\i <e) <^ SeAr^ = 2aie 
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We note that this estimate is almost exact, since the expansion factor is 
almost constant on each Wi due to (j5.4|) . 

Next, to obtain an open {60, l)-subset C as required in ()4.5|) . we 
need to subdivide the intervals Wl (Z into subintervals whose F- images 
are shorter than 60. More precisely, let us divide each curve F{Wl) whose 
length exceeds 60 into ki equal subintervals of length between 60/ 2 and 60, 
with ki < 2\F{Wl)\/5o. If \F{Wl)\ < So, then we set h = and leave 
unchanged. Then the union of the preimages of the e-neighborhoods of the 
new partition points on the curves F(Wl) has i^^y-measure bounded above 
by 

i i 

where we increased the numerical coefficient from 6 to 7 in order to take 
incorporate the factor e^^^°^ resulting from the distortion bounds (j5.4j) . Thus 
we obtain 

(5.7) ^w{rvl,i < 2aie + 7e6o^\W\ 

Lemma is proved. □ 

The one-step expansion estimate ()5.5|1 is geometrically explicit and easy- 
to-check, but unfortunately, as stated, it fails too often. Indeed, suppose 
5* C 5 is a singularity curve that divides an unstable manifold W into two 
components, Wi and W2- For ()5.5|) to hold, we need A]~^ + ^ < 1, which 
is a fairly stringent requirement (it means, in particular, that Ai > 2 or 
A2 > 2). What if the hyperbolicity of the system at hand is not that strong, 
i.e. the expansion of unstable curves is less than twofold? Then (15. 5|) fails, 
and so does (j4.5|) . 

In that case, according to the remark after Theorem it will be enough 
to prove ()4.5|) for any power F™. If we apply Lemma |H1 to the map F"^ we 
get the following: 

Lemma 9. The analogue of the condition for is equivalent to 

(5.8) am- =liminf sup 'S^A^^<1, 

where the supremum is taken over unstable manifolds W and Am,i, « > 1, 
denote the minimal local expansion factors of the connected components of 
W\Sm under F"^. 
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We call (|5.8|) the m-step expansion estimate. 

To summarize our discussion, we state another theorem: 

Theorem 10. Let F be defined on a 2D manifold M. and satisfy the basic re- 
quirements Suppose the singularity set S has the structure described 
by 115.1]} or by Remark S, and let i5.t^} hold. Then (a) there is a horseshoe 
Aq C M such that the return times R{x] F, Aq) satisfy the tail bound \3. 1]) . 
Thus, (b ) the map F : M ^ M enjoys exponential decay of correlations. 

In our applications, 4.1-4.6 will always hold and the singularity set S will 
obviously have the structure described by (jS.ip or by Remark S. The only 
condition we will need to verify is ()5.8|) for some m > 1. But still, it involves 
a higher iteration of F, which might be practically inconvenient, so we will 
simplify this condition further in the next section. 

6 A practical scheme 

The verification of (j5.8|) can be done according to a general scheme that we 
outline here. 

Suppose Sm is a finite union of smooth compact curves that are uniformly 
transversal to unstable manifolds (in our applications, the tangent vectors 
to those curves lie in stable cones ensuring transversality) . For an unstable 
manifold W, let Km{W) denote the number of connected components of 
W \ Sm . We call 

Km = liminf sup Km{W) 

So^O W: \W\<So 

the complexity of the map F™. Note that Km does not depend on the number 
of the curves in 5^, which may grow rapidly with m, but rather on the 
maximal number of those curves that meet at any one point of M. 

For example, on Fig. El a set of 10 curves of Sm is depicted; to ensure 
transversality, those curves run in a general vertical direction, while unstable 
cones (depicted by a shadowed triangle on the right) are assumed to be nearly 
horizontal; the complexity Km here equals 3, since sufficiently short unstable 
manifolds intersect at most 2 curves in Sm- 

Lemma 11. Suppose that for some m > 1 

(6.1) i^m<A"^ 
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Fig. 3: Singularity curves and unstable cones. 

where A > 1 is the minimum expansion factor of unstable vectors. Then 
/ TO) holds, thus (O holds for F"^. 

Proof. By the chain rule, we have j > A"* in ()5.8p . hence the lemma. □ 

The condition ()6.H1 and alike are known as complexity bounds in the 
literature. In some cases, the sequence Km has been proven to grow with 
m at most polynomially. Km = 0{m^) for some A > 0, thus (jfi.lll holds for 
all large enough m. It is believed that for typical billiard systems Km grows 
slowly (at least more slowly than any exponential function), but only in few 
cases exact proofs are available. 

The above arguments are valid if Sm is a finite union of smooth com- 
pact curves. However, in many billiards systems (and in all our applications) 
Sm consists of countably many smooth compact curves, which accumulate at 
some places in the space M. For example, the boundaries of homogeneity 
strips introduced in 4.4, converge to two lines ip = ±7r/2, hence their preim- 
ages (included in S) converge to the preimages of the lines ip = ±7r/2. Note 
that since Sm is closed, the accumulation points of singularity curves belong 
to Sm as well (they may be single points or some curves of Sm)- 

In this case the singularity curves S of the map F can be usually divided 
into two groups: a finite number of primary curves and finitely or countably 
many sequences of secondary curves, each sequence converges to a limit point 
on a primary curve or to a primary curve in S. In each example, the partition 
into primary and secondary singularity curves will be quite natural. 

Now let ly be a short unstable manifold that crosses no primary singular- 
ity curves but finitely or countably many secondary singularity curves. Then 
W\S consists of finitely or countably many components Wi, and we denote 
by Aj, as before, the minimum local expansion factor of Wi under the map 
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F. In all our applications, the accumulation of secondary singularity curves 
is accompanied by very strong expansion, so that Aj are very large. The 
following is our crucial assumption: 

Assumption (on secondary singularities). We have 

(6.2) =liminf sup VA^^ < 1, 

So^O W:\W\<So~' 

where the supremum is taken over unstable manifolds W intersecting no 
primary singularity curves and Aj, z > 1, denote the minimal local expansion 
factors of the connected components oiW \ S under the map F. 
We denote 

Oi: = max{^o,A"^} < 1. 

For the primary singularity curves, we assume an analogue of the com- 
plexity bound (|6.ip as follows. Let Sp be the union of the primary singularity 
curves and Spm = SpUF^^SpU- ■ -UF^^+^cSp. For an unstable manifold W, 
let Kp^miy^) denote the number of connected components of W \ Sp^m- i-e. 
the number of components into which W is divided by primary singularities 
only during the first m iterations of F. We call 

-^p,m = liminf sup Kp^miW) 

<5o^O \w\<So 

the complexity of the primary singularities of F"^. We now make our last 
assumption: 

Assumption (on primary singularities). For some m > 1 

(6.3) Kp,^ < 9^"^ 

Theorem 12. Assume \6. jjj) . // \6. 'J\i holds for some m > 1, then 45. 
follows, and thus \4.1^ is true for the map F™ . 

Proof. Let W CW he a. connected component of the set W \ Sp^m and Wj, 
j > 1, denote all the connected components of W \ Sm- Denote by A^- the 
minimum expansion factor of under the map F™. 

Lemma 13. In the above notation 

(6.4) liminf sup sup V[A'.]"i < ^7^. 

5o^O W:\W\<SoW'CW~^ 
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Proof. During the first m iterations of F, the images of W never cross any 
primary singularity curves but may be divided by secondary ones into finitely 
or countably many pieces. 

Now the proof goes by induction on m. For m = 1, the statement follows 
from ()6.2|) if W intersects S, or from the mere hyperbolicity 4.2 otherwise. 
Next we assume (j6.4j) for some m and apply this same argument to each 
connected component of the set F^{W') and then use the chain rule to 
derive ()6.4j) for m + 1. □ 

Lemma IT^ and the assumption ()6.3|) imply 

am < Kp^mOT < 1 

which proves the theorem. □ 

In the following sections, we apply our methods to three classes of chaotic 
billiards with slow mixing rates. 



7 Proof of Theorem 1 

Let 7?. C be a rectangle and Bi,...,Br C intTl open strictly convex 
subdomains with smooth or piecewise smooth boundaries, as described in 
Theorem 1. 

The billiard system in Q = TZ\UiBi is semidispersing. The collision space 
can be naturally divided into two parts: dispersing 

M+ = {(r,(^): r G UidB,} 

and neutral 

Mo = {(r,(^): r G dU} 

We set M = and consider the return map F : M ^ M. Thus we skip 
all the collisions with the flat boundary &TZ when constructing F. 

The map F can be reduced to a (proper) collision map corresponding 
to another billiard table by a standard "unfolding" procedure. Instead of 
reflecting a billiard trajectory at dTZ we reflect the rectangle TZ with all the 
scatterers Bi across the side which our billiard trajectory hits. Then the 
mirror image of the billiard trajectory will pass straight to the new copy of 
TZ (as if the boundary dTZ were transparent), see Fig. EJ 
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Fig. 4: Unfolding of a billiard trajectory. 

The mirror copies of TZ obtained by successive reflections across their sides 
will cover the entire plane M^. The new (unfolded) billiard trajectory will only 
hit the scatterers Bi and their images obtained by the mirror reflections. It is 
clear that those images make a periodic structure in M^, whose fundamental 
domain TZi consists of four adjacent copies of TZ, see Fig. HI In fact, TZi can be 
viewed as a 2-dimensional torus, with opposite boundaries identified. Thus 
we obtain a billiards in a table Qi C TZi with some internal convex scatterers 
and periodic boundary conditions. Our map F reduces to the proper collision 
map in this new billiard table. 

The billiard table Qi is dispersing since all the scatterers are strictly 
convex. It has "infinite horizon" (or "unbounded horizon"), because the free 
paths between successive collisions may be arbitrarily long. We call x G M 
an infinite horizon point (or, for brevity, IH-point) if its trajectory forms a 
closed geodesic on the torus TZi that only touches (but never crosses) the 
boundary dQi, and at all points of contact of this geodesic with dQi the 
corresponding components of dQi lie on one (and the same) side of this 
geodesic. There are at most finitely many IH-points x e M, and we denote 
them by Xq, g > 1, they play a crucial role in our analysis. 

First, let us assume that the scatterers Ei have smooth boundaries. Then 
the billiard in Qi is a classical Lorentz gas, also known as Sinai billiard table. 
The exponential decay of correlations (in fact, both statements (a) and (b) 
of Theorem for this billiard are proved in P|, Section 8. But the argument 
in PI is rather model-specific. Below we outline another argument based 
on our general scheme described in the previous sections. This outline will 
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also help us treat piece-wise smooth scatterers Bi later. We mention some 
standard properties of dispersing billiards and refer the reader to [71 |Hl E] 
for a detailed account. 

The discontinuity curves of the map F (the "old" singularities, in the 
terminology of Section 4.4) are of two types. There are finitely many long 
curves dividing M into finitely many large domains. In addition, there are 
infinite sequences of short singularity curves that converge to the IH-points 
Xq = {rq,ipq) G M (uotc that (fq = ±7r/2). The singularity curves divide 
the neighborhoods of the IH-points into countably many small regions (com- 
monly called cells). 

The structure of singularity curves and cells in the vicinity of every IH- 
point Xq is standard - it is shown on the Fig. El (a). There is one long curve 
Sq running from Xq into M and infinitely many almost parallel short curves 
Sq,n, n > 1, running between Sq and the border dM and converging to Xq 
as n goes to oo. The dimensions of the components of this structure are 
indicated on Fig. El (a). We denote by Mg„ the domain bounded by the 
curves Sq^n, Sq^n+i, Sq and dM (this domain is often called n-cell). 
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Fig. 5: Singularity curves and cells near an IH-point. 

We declare all the short singularity curves Sq^n with n > rii secondary, 
where ni is sufficiently large and will be specified below. All the other old 
singularity curves are declared primary. 

The "new" singularities (the preimages of the boundaries of the homo- 
geneity strips Hk, see Section 4.4) consist of infinite sequences of curves that 
converge to the old singularities. In particular, every curve Sq^n is a limit 
curve for an infinite sequence of "new" singularity curves, see Fig. El (b). All 
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the new singularity curves are declared secondary. 

It is known that the complexity of the primary singularities of the 
map F"^ grows at most linearly, i.e. Km < Ci + C2m, where Ci,C2 > 
are constants (see [B], Section 8). This means that no more than Ci + 
singularity curves of Sm meet at any one point x G M. This implies fj6.3p for 
any < 1 and sufficiently large m. 

It remains to verify the main assumption (j6.2|) . Suppose first that a short 
unstable manifold W crosses new singularity curves converging to a primary 
old curve. Then W is divided into countably many pieces Wk = WnF~^{Hk). 
The expansion of Wk under F is bounded below by > C/c^, where C > 
is a constant, see P, Equation (7.3). Thus 

oo oo 
ko ko 

which can be made < 1/2 by choosing ko large enough, say ko > AC~^. 

Next, let a short unstable manifold W intersect some (or all) secondary 
old curves Sq^n with some q > 1 and rii < n < oo, and near each of them 
it intersects infinitely many new singularity lines converging to the old one. 
For each n > rii denote by Wn the piece of W between Sq^n and Sq^n+i, 
and let Wn,k = Wn H F~^{Hk) for n > rii and k > ko denote the connected 
components ofW\S. It is important to observe that the image F{Wn) only 
intersects homogeneity strips with k > x^^^-, where x > is a constant 
(see in], page 544). The expansion of Wn,k under F is bounded below by 
An,fc > Cnk? ^ where C > is a constant, see [0], Equation (8.2). Thus we 
have 

oo oo oo oo 

E E <- E E 



n=ni 



11-1/4 



<ioc-VX 

which can be made < 1/2 by choosing rii large enough, say ni > [20 C^^x'^T- 
This concludes the verification of ()6.2|) . By Theorem^] and the remark 
after Theorem the assumptions ()4.5|) - ()4.7|) hold, hence the induced map 
F: M ^ M has exponential mixing rates by Theorem El 
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Lastly, in order to determine the rates of mixing for the original collision 
map J-": Ai ^ Ai of the billiard in Q = 7^ \ UjSj, we need to estimate the 
return times defined by ()3.6|) . If the short singularity curves Sq^n are labelled 
in their natural order (so that n increases as they approach the limit IH-point 
Xq), then 

(7. 1) {x e M: R{x;J',M) > n} CUq Ui>,„ M,,, 

where c > is a constant and Mq^i denotes the z-cell. This is a simple 
geometric fact, we leave the verification to the reader. 

It is easy to see on Fig. El (a) that /i(Mq^„) = 0{n~^), because the 
"width" of Mq^n (its r-dimension) is 0{n~'^), its "height" (the yp-dimension) 
is 0(n~^/^) and the density of the /i measure in the (r, coordinates is 
cos(p = 0(n~^/^). Hence /i(i?(x; JF, M) > ra) = 0{n~'^). Theorem |3] now 
implies the required bound on correlations 

|Cn(/,fl',^,yu)| < const ■ {Innf/n. 

To complete the proof of Theorem 1, we need to consider scatterers Bi 
with piece-wise smooth boundaries. Now there are two types of IH-points: 
those whose trajectories intersect (graze) dQi at points where dQi is smooth 
we call them IHl-points, and those whose trajectories intersects dQi at corner 
points, as shown on Fig. we call them IH2-points. 

The analysis of IHl-points is identical to the one above. The structure of 
singularity curves and cells in the neighborhood of an IH2-point Xq = [vq, ipq) 
is shown on Fig. ^ (b), we refer the reader to [B], Section 4, for more detail. 
It is important to note that \ipq\ < 7r/2 and the first few F-images of a small 
neighborhood of Xp do not intersect dM (they "wander" inside M, as it is 
explained in jH]), so that there is no "new" singularity curves inside the cells 
Mq^n for large enough n. 

Following [7], we assume that the condition ()6.3p holds for some m > 1 
and proceed to verify the main assumption ()6.2p . Let /Cmin > denote the 
minimum curvature of the boundary of the scatters. Let d denote the width 
(the smaller dimension) of TZ and let 

b = mm(dist{Bi, Bj), dist{Bi,dn)). 

Let ly be a short unstable manifold and denote by Wn the piece of W between 
Sq^n and Sq^n+i- Unstable manifolds are represented by increasing curves in 
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Fig. 6: Discontinuity curves and cells near an IH2-point Xg. 



the rip coordinates. It is easy to see on Fig. IHl (b) that unstable manifolds 
cannot cross infinitely many singularity curves. In fact, by a simple geometric 
calculation, there exists a constant 



1 < a < 1 



^min ~l~ 



cos <pq 

2b 



k 



mill 



cos ipq 



such that if W crosses the curves Sq^n for all n ^ [ni,n2\, then n2 < arii. 
The expansion of Wn under F is bounded from below by 



An > 1 



Cn, 



where C = JCmmd/ cos(^q, see [0], Equation (6.8). Here we used the fact that 
the intercollision time for points x G Wn is at least dn. Thus we have 



an ^ ^ an ^ 

y—<-y- 



(7.2) 



< 



C + 7] 



where r) = d/2b and we used the fact that ln(l + x) < x for a; > 0. 

The bound ()7.2|) is finite, but we need it to be less than one, which only 
happens if ?7 < C^, that is if cos^ ipg < 2bdlC'^^^. While this is indeed the case 
for some billiard tables, it is not hard to enforce the condition ()6.2|) for all 
relevant billiard tables by considering a higher iteration of the map F, as it 
is explained in [7j, page 104. 
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By Theorem ^1 and the remark after Theorem |3 the assumptions (|4.5p - 
(|4.7|) hold, hence the induced map F : M ^ M has exponential mixing rates. 

Lastly, in order to determine the rates of mixing for the original collision 
map JF: — > 7W of the billiard in Q = 7^ \ Uji3j, we need to estimate the 
return times defined by ()3.6|) . Similarly to ()7.1|) . we now have 

{xeM: R{x; M) > n} C U, \^m>c^n M^^^n 

where Ci > is a constant. 

It is easy to see on Fig.|ni(b) that fi{Mq^n) = 0{n~^) for IH2-points, hence 

ld{R{x;J^,M) >n) = 0{n-^). 

Theorem 13 now implies the required bound on correlations 

|C„(/,fi', JF,yu)| < const ■ (\nnY/n. □ 

8 Proof of Theorem 2 

Let Q C be a billiard table satisfying the assumptions of Theorem 2. 
Its boundary can be decomposed as dQ = d^Q U d~Q so that each smooth 
component Fj C d^Q is dispersing and each component Tj C d^Q is focusing 
(an arc of a circle that is contained in Q). The collision space can be naturally 
divided into two parts: 

M± = {{r,ip) : r G d^Q}. 

As distinct from dispersing billiards, expansion and contraction in A4_ 
are not uniform. More precisely, the expansion and contraction are weak (per 
collision) during long sequences of successive reflections in the same focussing 
boundary component. 

These sequences become a "disturbing factor" similar to reflections in 
neutral (fiat) boundary components studied in the proof of Theorem 1, hence 
the induced map F: M — > M must be defined so that those sequences are 
skipped. We set 

(8.1) M = M+U{x e M-: 7i{x) e r„ 7r{J^-^x) e Tj, j ^ i}, 

where 7t{x) = r denotes the first coordinate of the point x = (r, ip). Observe 
that M includes all collisions at dispersing boundary components but only 
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the first collision at every focusing component. We note that in [71 the 
last collision (rather than the first) in every arc is used in the construction of 
an induced map, but that definition of M leads to unpleasant complications 
in the analysis, which we will avoid here. 

There are two ways in which billiard trajectories can experience arbitrarily 
many reflections inside one focusing component (arc) Fj C d~Q. First, if the 
collisions are nearly grazing [ip ^ ±7r/2), the points x, T{x), jF^(x), . . . are 
close to each other, so that the sequence {jF"(x)} moves slowly along the arc 
Fj until it comes to an end of the arc and escapes. This is possible for arcs 
of any size. 

Second, if y9 ^ 0, then the billiard trajectory runs near a diameter of 
Fj, hits Fj on the opposite side and then comes back, so that the points x, 
J-''^{x), J-'^{x), . . . are close to each other. Then the two sequences {^^""(x)} 
and {J-''^"'^^{x)} move slowly along the arc Fj until one of them finds an 
opening in Fj and escapes. Similarly, the trajectory can run close to a periodic 
trajectory inside Fj of any period p > 2. All these, however, require the arc 
Fj to be larger than half-circle, which is specifically excluded by assumption 
(ii) of Theorem 2. So we do not need to deal with these trajectories now, 
but we will encounter them later. 

In the coordinates (r, (p), the set M is the union of rectangles and cylinders 
corresponding to the dispersing boundary components and parallelograms 
corresponding to the focusing boundary components, as the one shown on 
Fig. [7| In each parallelogram, unstable manifolds are represented by decreas- 
ing curves in the (r, ip) coordinates, and singularity manifolds - by increasing 
curves (which is opposite to dispersing billiards). Between successive colli- 
sions at the focusing boundary, each unstable manifold first converges (fo- 
cuses, or collapses), passes through a conjugate (defocusing) point, and then 
diverges. We refer the reader to [SI 111^3 for a detailed account of hyperbol- 
icity in Bunimovich billiards. It is important that the conjugate point always 
lies in the first half of the segment between the consecutive reflections, which 
guarantees a monotonic expansion of unstable manifolds from collision to 
collision in a special p-metric on unstable curves, defined by dp = cos (f dr, 
see ini |7j for definitions and details. 

The induced map F has infinitely many discontinuity lines, which accu- 
mulate in the neighborhood of two vertices (top and bottom) of each paral- 
lelogram, see Fig. [7| Each such vertex Zq is a limit point for infinitely many 
almost parallel straight segment Sg^n C S running between two adjacent 
sides of the parallelogram. They divide the neighborhood of Zg into count- 
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Fig. 7: Sliding trajectories and cell structure of Bunimovich billiard tables 

ably many cells, which we denote, as before, by M^^, n > 1. Here the cell 
Mq^n consists of points experiencing exactly n consecutive collisions with the 
arc (counting the first collision included in M). In fact, for n large enough, 
those trajectories are more like "sliding" along the arc rather than reflecting 
off it, hence the corresponding singular points are said to be of sliding type, 
see more details in jH] and [7j. 

We declare all the short singularity curves Sg^n with n > rii secondary, 
where rii is sufficiently large and will be specified below. All the other sin- 
gularity curves are declared primary. 

Following [Zj, we assume that the condition holds for some m > 1. 
We remark that ()6.3|) holds for generic billiard tables, i.e. for an open dense 
set of billiard tables in the metric, which one can show by using standard 
perturbation techniques, but we do not pursue this goal here. 

We now proceed to verify the main assumption (j6.2|) . Let W C Mg „ be 
an unstable manifold, so that F = JF" on W. Then the map JF"~^ expands 

by a factor 0{n), and then JF expands the manifold J-'"'~^{W) by a factor 
0{n), see jTj, Section 2.5. Thus, the induced map F expands W hj a. factor 
> CTi^, where c > is a constant, and we have 



which can be made < 1/2 by choosing ni large enough, say rii > 4/c. 

By Theorem El and the remark after Theorem the assumptions ()4.5|1 - 
()4.7p hold, hence the induced map F : M ^ M has exponential mixing rates. 

Lastly, we need to estimate the return times (3.6) to determine the rates 
of mixing for the original collision map J^: M. ^ M.. The cell M^ ^ has 
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"width" "height" see Fig.[71 and the density of the measure 

/i on Mg^n is 0{n~^), hence fi{Mg^n) = Since 

{xeM: R{x; M) > n\ C U^>„ M<j,™, 

we get ii{R{x\T^M) > n) = (9(n~^). Theorem 4 now imphes the required 
bound on correlations 

|Cn(/,fi',-^,yu)| < const ■ {Vnnflr?. 

Note that the faster rate of the decay of correlations here (as compared to 
Theorem 1) is purely due to the smaller cells Mg^„, and not due to any 
properties of the induced map F . □ 



9 Proof of Theorem 3 

Let Q C be a stadium. Its boundary can be decomposed as 

dQ = d°Q U d-Q, 

where dQ^ = riUr2 is the union of two straight sides of Q, and dQ~ = r3Ur4 
is the union of two arcs. The collision space can be naturally divided into 
focusing and neutral parts: 

Mo = {{r,ip) : r e d^Q}, M. = {{r,ip) : r E O-Q}. 

We define the induced map F: M — > M on the set 

(9.1) M = {xeM^: 7r{x) G Ti, n{J^~^x) G Tj, j ^ i}, 

Note that M only contains the first collisions with the arcs (the collisions 
with the straight lines are skipped altogether). In the coordinates {r,if) the 
set M is the union of two parallelograms. 

First let us consider the straight stadium, see Fig. 121(c). This model was 
already handled by Markarian ^Hj, but we present here a simplified proof 
based on our general method. 

The map F can be viewed as a collision map corresponding to another 
billiard table obtained by "unfolding" the stadium along the straight bound- 
ary segments. In other words, instead of reflecting a billiard trajectory at 
d^Q, we reflect the stadium across its straight boundaries, see Fig. |H1 
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Fig. 8: Unfolding of straight stadium 



Then the mirror image of the bilhard trajectory will pass straight through 
dQ until it meets the mirror image of the arcs and get reflected there. The 
mirror copies of the straight stadium obtained by successive reflections about 
their straight sides make an unbounded strip, see Fig. |H1 The new billiard 
in the unbounded strip has "infinite horizon", since the free path between 
collisions may be arbitrarily long. 

The singularity set iS C M of the map F consists of two types of infinite 
sequences of singularity curves, as shown on Fig. |H1 The first type accumulate 
near the top and bottom vertices of the parallelograms, they are generated 
by trajectories nearly "sliding" along the circular arcs. These are identical 
to the ones discuss in the proof of Theorem 2, so we omit them. 

The curves of the second type accumulate near the other two vertices 
of the parallelograms (those lie on the line </:> = 0), they are generated by 
trajectories experiencing arbitrary many bounces between the two straight 
sides of the stadium (i.e. by long collision-free flights in the unbounded strip). 
The structure of those singularity curves is shown on Fig. El (b), see 0, 
Section 6.3, for more details. 

Actually, there are two types of cells near the vertices xi and X2 in M. 
The first type contains points mapped (by JF) directly to a straight side of 
the stadium, we denote them by M^, n > 1; the second type contains points 
that are mapped by T onto the same arc first and then to a straight side, we 
denote them by M^, n > 1. In Fig. |H1 the point z belongs to M| and F[z) 
belongs to M^. Points in and experience exactly n reflections off the 
straight sides before landing on the opposite arc of dQ. 

To prove the condition (j6.3j) it is enough to establish the linear growth of 
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yi 




Fig. 9: Discontinuity curves of straight stadium 



the complexity 

(9.2) K^^m <Ci + Cam, 

which is similar to the one obtained for the Lorentz gas in [H], Section 8 (recall 
our proof of Theorem 1). The proof of ()9.2j) for the Lorentz gas is based on 
the continuity of the billiard flow (as it is explained in ^U], Lemma 5.2). The 
flow in the stadium is obviously continuous as well, thus the same argument 
implies the linear estimate ()9.2|1 for the stadium, we omit details. 

We now proceed to verify the main assumption ()6.2j) . Our calculations 
are based on two known facts mentioned in [7j, Section 2.6, which can be 
verified by direct calculations. First, if an unstable manifold W intersects 
cells (and M^) with rii < n < n2, then 

n2 < 9ni + const. 

Second, the expansion factor of F on each piece W fl M^, n > 1, is > 4n and 
the expansion factor on each piece W fl is > 8n. Thus we have 

"2-1 "2-1 -| Q 

y -L< y f— + — )<-in9<i. 

Am ^ V4m 8m/ 8 

m=n\ m=ni 

By Theorem IT^ and the remark after Theorem the assumptions ()4.5p - ()4.7p 
hold, hence the induced map F: M ^ M has exponential mixing rates. 

Next we need to estimate the return times (3.6) to determine the rates 
of mixing for the original collision map J^: M. ^ Ai. The n-cells near the 
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top and bottom vertices yi and ?/2 have measure (9(n~^), see the previous 
section. The n-cells and have measure 0{n~^), as one can easily 
see from Fig. IHl (b) (note that ~ 0, hence cos 93 ~ 1). Thus we get 
li(^R{x; J-', A4) > n) = 0{n~'^). Theorem 4 now imphes the required bound 
on correlations 

< const ■ {Inny/n. 

Note that there are two types of cells here, and the rate of the decay of 
correlations is determined by the "worse" (larger) cells and (those 
near the points xi and X2). 

This concludes the proof of Theorem 3 for the straight stadium. □ 

Next we consider a skewed stadium (or a "drive-belt" table), see Fig. ITfll 
(a). Again, we can unfold the skewed stadium by reflecting it repeatedly 
along the fiat boundaries, as shown on Fig. El (b). The new billiard table 
has similar structure to the Bunimovich type billiard tables of the previ- 
ous section, but it does not satisfies assumption (ii) of Theorem 2, since it 
necessarily contains an arc larger than half a circle. 




(a) (6) 

Fig. 10: Unfolding of a skewed stadium 

The "unfolding" of skewed stadium also shows that it has "finite horizon" , 
since the free path between successive collisions is uniformly bounded from 
above, thus there are no cells like and described above, but there is 
a new, equally influential, type of cells, see below. 
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In the coordinates (r, ip) the set M is the union of two parallelograms, 
one corresponds to the smaller arc, the other to the larger arc, we call them 
small and hig parallelograms, respectively. 

The singularity set 5 C M of the map F consists of two types of infinite 
sequences of singularity curves. Curves of the first type accumulate near 
the top and bottom vertices of each parallelogram, they are generated by 
trajectories nearly "sliding" along the circular arcs. These are identical to 
the ones discuss in the proof of Theorem 2, so we omit them again. 

Curves of the second type accumulate near two points xi and X2 on the 
slanted sides of the big parallelogram (where those sides intersect the line 
^9 = 0). There are two parallel long singularity lines Si and 5*2 starting at 
Xi and X2, respectively, and running into M, and two infinite sequences of 
almost parallel straight segments Sq^n, where g = 1,2 and n > 1, running 
between Sq and the nearby side of M and converging to as — > oo, see 
Fig. EH 

The singularities of this second type are generated by trajectories expe- 
riencing arbitrary many collisions with the large arc while running almost 
along its diameter (this type of trajectories was described in the previous 
section, but not studied in detail there, because it was ruled out by the as- 
sumption (ii) of Theorem 2). Denote by Mg^„ the n-cell bounded by Sq^m 
Sq^n+i, Sq aud dM. The n-cell Mq^n consists of points experiencing exactly 
n collisions with the large arc. 



c/n 














Si 




M 



Fig. 11: The cell structure of skewed stadium 

There is a peculiar feature of the new cells Mg „. Points in the upper half 
of Mq^n near Xq (the dark grey area marked by y on Fig. ITTj) are mapped 
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by F into the cells M3_gm; "^i < m < n2, near the other limit point X3_q, 
and n2 < 49?ii + const (see the dark grey area marked by Fy). However, 
points in the lower half of the cell Mq^n (the light grey area marked by z on 
Fig. [TT|) are mapped under F into a neighborhood of the other limit point 
Xs-q above or below the line Ss-q without crossing any singularity lines there 
(see the light grey area marked by Fz), and the second iterate of F, as it is 
illustrated by Fig. (b), maps them into the cells M3_q m, rii < m < n2 
with n2 < 49ni + const (see the light grey area marked by F^z). The bound 
712 < 49ni + const follows from elementary (but tedious) calculations which 
we leave out. 

Fig. ^1 (b) shows the dynamics of the lower half of the cell Mg „. Any 
point z in that part experiences n collisions with the large arc, then crosses 
a straight side of Q and lands on the adjacent copy of the large arc (this 
becomes Fz), but then it crosses the same straight side back and lands on 
the old large arc again, where it starts another long series of m reflections. 
Hi < m < 71.2, running nearly along the arc's diameter. 

The above analysis implies that if an unstable manifold W intersects cells 
Mq„ with Hi < 772, then n2 < 49ni + const. Next, by a direct calculation 
(we omit details) the expansion factor of the map F on the curve W (1 Mq^n 
is > 8n. However, since Wn = W H Mq^n consists of two parts (the upper 
half, call it W^, and the lower half, call it W^), which evolve differently, the 
number of pieces of W in our estimate is doubled and we get 



This finishes the proof of ()fj.2|l . Now by Theorem El and the remark af- 
ter Theorem the assumptions ()4.5p - ()4.7p hold, hence the induced map 
F : M — s> M has exponential mixing rates. 

Lastly, we need to estimate the return times (3.6) to determine the rates 
of mixing for the original collision map !F : A4 ^ Ai. The cells near the 
top and bottom vertices of the parallelograms of M have measure 0{n~'^), 
see the previous section. The n-cells Mn have measure 0{n~^), as one can 
easily see from Fig. ^2 (note again that ip ^ 0, hence cosip ^1). Thus we get 
^[R{x]T,M.) > nj = 0{n~^). Theorem 4 now implies the required bound 
on correlations 




- In 49 < 1. 

4 



Cn{f,g,^,fJ')\ < const ■ (Innf/n. 



The proof of Theorem 3 is complete. 



□ 
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10 Other examples and open questions 



Here we discuss two types of chaotic billiard tables for which our method 
fails. This discussion will demonstrate the limitations of the method, in its 
present form, and indicate directions for future work. 

First, let Q be a Bunimovich type table satisfying the assumption (i), but 
not (ii), of Theorem 2, that is let d~Q contain an arc Fj larger than half a 
circle. Suppose we define the induced map T: M ^ M by 1)8.11) again, then 
M will contains a parallelogram Mi corresponding to the arc Fj, and the 
structure of the singularity lines in Mj will be the same as the one shown on 
Fig. ^2 In particular, the expansion factor of F on any unstable manifold 
W n Mq „ will be ~ cn, where c > is a constant. 

However, unlike the case of a drive-belt region illustrated on Fig. ITU 
now some unstable manifolds may intersect infinitely many cells Mg^„, in 
particular all n-cells with n > rii for some ni > 0. In that case 



and so the condition ()fi.2j) fails. And it fails in a major way, since neither F 
nor any power F"^ can possibly satisfy ()6.2|) . 

This failure poses interesting questions. Does this mean that the induced 
map F: M — > M has subexponential decay of correlations? If not, can 
the method be improved to overcome the trouble and establish exponential 
mixing for F? Or should one seek a different definition of the induced map F 
in order to avoid the trouble? These questions remain open at the moment. 

A similar failure occurs for a special modification of the stadium, where Q 
is bounded by two parallel straight segments and two arcs which are less than 
half a circle (in that case the arcs have to be transversal to the segments at 
the intersection points). For the hyperbolicity and ergodicity of this model, 
see 0, Section 6.3. Suppose we define the induced map F: M — > M by 
(jnnj)) as before. Then one can investigate the singularities of this map and 
find (see Section 6.3 of for more details) that an unstable manifold W can 
be broken by S into infinitely many pieces Wn, n > rii, and the expansion 
factor of the map F on Wn is ~ cn, where c > is a constant. Hence again 
we arrive at a divergent series similar to (|lU.ip . so the condition (|6.2p fails. 
It is not clear what this trouble implies for the map F and how to deal with 
it. 



(10.1) 




36 



These two example demonstrate the hmitations of our algorithm, in its 
present form. There are other models, not discussed in this paper, to which 
our scheme probably applies but it requires substantial extra effort. One of 
them is a dispersing billiard table with cusps (corner points where the ad- 
jacent boundary components are tangent to each other). The estimation of 
mixing rates for this model remains an open problem since the first publi- 
cation ^7] on it in 1983. The other is a dispersing billiard table where the 
curvature of the boundary vanishes at some points, i.e. where the boundary 
looks like the graph of function y = c\x\^ near x = with some (3 > 2 (then 
the curvature vanishes at x = 0). Our preliminary calculations show that 
the rate of mixing is polynomial and its degree depends on j3. The work on 
both models is currently underway. 

Another interesting question is how to bound the correlations from below 
(recall that all our bounds are from above). It is intuitively clear from the 
results of that the polynomial bounds on correlations we established here 
(up to the logarithmic factor, though) are sharp, i.e. cannot be improved. 
However, formal lower bounds on correlations may be hard to obtain. We 
only point out recent results by O. Sarig in this direction pU] . 
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